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ABSTRACT 

The hard X-rays in Cyg X-l and similar black hole sources are possibly produced in an active corona atop an 
accretion disk. We suggest that the observed weakness of X-ray reflection from the disk is due to bulk motion 
of the emitting hot plasma away from the reflector. A mildly relativistic motion causes aberration reducing X- 
ray emission towards the disk. This in turn reduces the reprocessed radiation from the disk and leads to a hard 
spectrum of the X-ray source. The resulting spectral index is V w 1.9/vB where B = 7(1 + j3) is the aberration 
factor for a bulk velocity (3 = v/c. The observed V w 1.6 and the amount of reflection, R w 0.3, in Cyg X-l 
in the hard state can both be explained assuming a bulk velocity (3 ~ 0.3. We discuss one possible scenario: the 
compact magnetic flares are dominated by e ± pairs which are ejected away from the reflector by the pressure of 
the reflected radiation. We also discuss physical constraints on the disk-corona model and argue that the magnetic 
flares are related to magneto-rotational instabilities in the accretion disk. 

Subject headings: accretion, accretion disks - black hole physics - gamma-rays: theory - radiation mechanisms: 
thermal - stars: individual: Cyg X-l - X-rays: general 



1. INTRODUCTION 

Galactic black holes (GBHs) in their hard state and radio- 
quiet active galactic nuclei (AGNs) have similar X-ray spec- 
tra which are well explained by Comptonization of seed soft 
photons in a hot cloud with scattering optical depth tt ~ 1 
and temperature kT ~ 100 keV (see reviews by Zdziarski et 
al. 1997; Poutanen 1998). The Compton reflection feature in 
observed spectra and a fluorescent iron line indicate the pres- 
ence of relatively cold gas in the vicinity of the X-ray source. 
The most likely reflector is an accretion disk, and the hard X- 
rays are possibly produced in a hot corona of the disk (e.g., 
Bisnovatyi-Kogan & Blinnikov 1977; Liang 1979). The dis- 
sipation of magnetic energy in the corona may feed the X-ray 
luminosity (Galeev, Rosner, & Vaiana 1979), while the underly- 
ing accretion disk reprocesses the incident X-rays and supplies 
seed soft photons to the corona. 

A self-consistent disk-corona model was developed, in which 
the coronal plasma is cooled by its own radiation reprocessed 
in the disk (Haardt & Maraschi 1993). The slope of the emerg- 
ing X-ray spectrum is then determined by the fraction of the 
X-ray luminosity which comes back to the source as a repro- 
cessed soft radiation. In a self-consistent model, this fraction 
equals A^ 1 where A is the Compton amplification factor of the 
hot corona. The resulting spectrum depends on the geometry 
of the corona. Ovservations favour a patchy corona consisting 
of separate active blobs (Haardt, Maraschi, & Ghisellini 1994; 
Stern et al. 1995; Poutanen & Svensson 1996). 

The disk-corona model has recently been found to be in 
conflict with observations of Cyg X-l and similar black hole 
sources in the hard state (e.g., Gierlihski et al. 1997): i) The ob- 
served spectrum is very hard which corresponds to a Compton 
amplification factor A <; 10 and implies soft photon starvation 
of the hot plasma (Zdziarski, Coppi, & Lamb 1990). The model 
predicts A ^ 5 unless the active blobs are elevated above the 
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disk at heights larger than the blob size (Svensson 1996). ii) 
The model with elevated blobs would yield a strong reflection 
component, R = fi/27r « 1, where SI is the solid angle cov- 
ered by the cold matter as viewed from the X-ray source. The 
reported amount of reflection is small, R ~ 0.3. 

The photon starvation and weak reflection may be explained 
if the cold reflector is disrupted near the black hole. This would 
agree with the idea that accretion proceeds as a hot flow in 
the inner region (Shapiro, Lightman, & Eardley 1976; Esin et 
al. 1998; Krolik 1998; Zdziarski 1998). One may fit the ob- 
served spectra with a toy model of inner hot cloud upscattering 
soft photons supplied by the surrounding cold disk or by dense 
cloudlets inside the hot region (e.g., Poutanen, Krolik, & Ryde 
1997; Zdziarski et al. 1998). 

However, the weak reflection does not necessarily imply that 
the inner cold disk is disrupted. One suggested alternative is 
that the apparent weakness of the reflection features may be 
due to a high ionization of the upper layers of the disk (Ross, 
Fabian, & Young 1998). A detailed analysis remains to be done 
to check whether the highly ionized reflector is consistent with 
observed spectra. 

In this Letter, we suggest another alternative. We find that 
the disk-corona model may be reconciled with observations if 
the plasma in the active regions has a mildly relativistic bulk 
velocity directed away from the disk. We argue that the plasma 
heated in magnetic flares is likely to be ejected with such a ve- 
locity, especially if it is comprised of e ± pairs. 

2. MAGNETIC FLARES IN THE CORONA 

The usually exploited model for the corona formation is 
that of Galeev et al. (1979, hereafter GRV). According to the 
model, a seed magnetic field is exponentially amplified in the 
disk due to a combination of the differential Keplerian rota- 
tion and the turbulent convective motions. The amplification 
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time-scale at a radius r is given by tc ~ r/3v where v is the 
convective velocity. GRV showed that inside luminous disks 
the field is not able to dissipate at the rate of amplification, and 
buoyant magnetic loops must emerge above the disk surface 
where the magnetic field may annihilate quickly. The rate of 
magnetic energy production per unit area of the disk equals 
Fb = 2hwB/tc where wb = B 2 /8ir is the magnetic en- 
ergy density and h is the half-thickness of the disk. Assum- 
ing that the magnetic stress t rip = B v B r /4ir is responsible for 
the transfer of angular momentum in the disk, one can com- 
pare Fb to the total surface dissipation rate, F + = 3t rip c s , 
where c s is the sound speed in the disk (Shakura & Sunyaev 
1973). One then finds F B /F + = h/r (taking into account 
that 2wB/t rip — B v /B r = c s /v in the GRV model). Hence, 
the GRV mechanism is able to dissipate only a small fraction 
<~ h/r <C 1 of the total energy. This is in conflict with the 
spectra of GBHs in the hard state, in which a large fraction of 
the energy is emitted in hard X-rays. 

Recent simulations of magnetic turbulence in accretion disks 
indicate that the magneto-rotational instability is likely to drive 
the turbulence (see Balbus & Hawley 1998 for a review). This 
instability operates on a Keplerian time-scale which is typically 
h/r times shorter than tQ. It produces magnetic energy with 
rate Fb ~ F+, i.e., fast enough to explain the bulk of energy 
release as magnetic dissipation. Combined with the GRV argu- 
ment for magnetic buoyancy, it follows that a large fraction of 
F+ may dissipate in the corona. 

The corona thus may be the place where magnetic stress driv- 
ing the accretion is transported to and released. For a standard 
radiation pressure dominated disk, t rip w m p cClK/^T, where 
fix is the Keplerian angular velocity. The magnetic field in the 
corona is likely to be strongly inhomogeneous and it is plausi- 
ble that the main dissipation occurs in localized blobs where the 
magnetic energy density, wb, is much larger than t rip given by 
the standard model. The accumulated magnetic stress may sud- 
denly be released on a time-scale to ~ 10r&/c (the "discharge" 
time-scale, see Haardt et al. 1994), where r& is the blob size. 
This produces a flare of luminosity L <~ r^ws/to, which may 
have a compactness parameter I = LaT/ri,m e c 3 ~ 10 — 10 3 . 

The compact flares can get dominated by e ± pairs created 
in 7 — 7 reaction (e.g., Svensson 1986). Pairs would keep an 
optical depth tt ~ 1 during the flare and upscatter soft radia- 
tion entering the e ± blob. The blob temperature is typically in 
the range 50 — 200 keV depending on I and on whether there 
is a non-thermal e ± tail. A strongly localized e ± flare with 
I ~ 10 2 — 10 3 may have a low temperature, kT < 100 keV, 
observed in some GBHs and AGNs. In the presence of a non- 
thermal tail, the flare may have even kT ^ 50 keV which is 
observed in the case of GX 339-4 (Zdziarski et al. 1998). 

Consider now one possible mechanism of bulk acceleration 
which should operate in an blob located above a reflec- 
tor. The blob luminosity is partly reflected, and hence the pair 
plasma is immersed in an anisotropic radiation field. The net 
radiation flux, F <~ L/r\, must efficiently accelerate the light 
pairs. The time-scale for acceleration to relativistic velocities 
is t a r~j m e c/f ~ l~ x r\,/c where / ~ Fcr^/c is the acceler- 
ating radiative force. The shortness of t a <C r^/c implies that 
the bulk velocity saturates at some equilibrium value limited by 
the radiation drag (e.g., Gurevich & Rumyantsev 1965; Sikora 
& Wilson 1981). We discuss the pair acceleration by reflected 
radiation in more detail in §4. 

The bulk motion in the hot blob does not necessarily mean 
that the blob itself moves. Instead, the dissipation region can 



be static. Then the escaping pairs are replaced by newly cre- 
ated e . The escaping pairs are immediately cooled down to 
the Compton temperature kTc ~ 10 keV. Note that in a bright 
7-ray source, 7 — 7 interactions would produce a lot of pairs 
above/between the active regions. An outflow of cold pairs cov- 
ering the whole inner region of the disk may then be created 
(Beloborodov 1998a). 

3. BULK MOTION AND REFLECTION 

It was noted previously that bulk motion of the flaring 
plasma may reduce the reflected radiation from an accretion 
disk (Wozniak et al. 1998). Reynolds & Fabian (1997) calcu- 
lated the impact of motion on the Fe Ka line. Here, we evaluate 
the amount of reflection and study how fast motion affects the 
coupling between the X-ray source and the reflector. We as- 
sume that the active blob luminosity is due to Compton ampli- 
fication of the reflected (reprocessed) radiation, and the multi- 
ply upscattered photons get isotropized in the plasma comoving 
frame. The angular distribution of the blob luminosity in the lab 
frame is then given by 



L(p) = 



(1) 



2 7 4(1-/3 M )3' 

where yu = cos 9, 9 is the angle between the ray and the plasma 
velocity, 7 = (1 - /3 2 )~ 1/2 , and L = j\ L(fj,)dfj, is the total 
luminosity of the blob. For a typical spectral index, T ~ 2, the 
relativistic transformation of the specific luminosity, L v (n), is 
the same as that of the bolometric luminosity (Rybicki & Light- 
man 1979). In equation (1), we assume that the dissipation re- 
gion is static and there is no retardation effect (cf. Rybicki & 
Lightman 1979). For a static blob with a life-time to ^ r b/c 
one may consider a time-independent picture in a fixed geome- 
try like the usual disk-corona model. The only difference is that 
now the hot plasma in the blob has a bulk velocity. 

We will assume that the velocity is perpendicular to the disk. 
The X-ray luminosity striking the disk is then given by 



,0 

L_ = / L(ji)dfi 



L(l + f3/2) 
2 7 4(l+/3)2 



(2) 



For a system inclination n 
reflection, R = Q./2-k, is 

R = 



cos 9, the apparent strength of 



(l+/3/2)(l-/3 M ) 3 
(1 + /3) 2 



(3) 



The reprocessed radiation from the disk is partly intercepted 
by the active blob. The intercepted soft luminosity may be ex- 
pressed in terms of an effective fi s , 



"-Ms 

X / L(/i)d/z. 



(4) 



The case /i s = describes a slab geometry of the active region, 
while ii s <~ 1/2 roughly corresponds to a blob with radius of or- 
der its height. The factor \ = 1 — a represents the efficiency of 
reprocessing, a w 0.1 — 0.2 is the disk albedo (e.g., Magdziarz 
& Zdziarski 1995). 

The fraction of L which returns as soft radiation is a function 
of (3 and fi s , 



L, = L sQ l + (l + /i,)/3/2 
L ~ L 7 4 (l + /3) 2 (l + ^/3) 2 ' 



(5) 
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where L s0 = \(1 — fi s )L/2 corresponds to the usual case of a 
static corona (0 = 0). With increasing 0, the intercepted soft 
luminosity is suppressed. 

The spectral index of the Comptonized radiation emerging 
from the blob is determined by the Compton amplification 
factor taken in the plasma comoving frame. A luminosity 
dL = L(p)dp transforms into the comoving frame as dL c = 
7 2 (1 — 0fi)dL = (1 + 0p, c )~ 1 dL, where index c stands for the 
comoving frame (see Rybicki & Lightman 1979). This gives 
L c = L and L c s = 7 2 (L S - 0® s ), where $ s = / L s (fi)fidfi 
is the total flux of the soft radiation through the blob. We ap- 
proximate L s (p) = L s j (1 — fj, s ) at (i > fi s and L s (/z) = at 
fi < n s . Then 



Then $diss = 0, and the equilibrium condition (8) combined 
with (6) yields the relation 



L s (l + fi s ) 



(6) 



which yields L C JL S = 7 2 [1 - 0(1 + £t s )/2]. We then get the 
amplification factor 



_L_ 2 7 2 (l + /3) 2 (l + Ms/ 3) 2 
L% x (l-^)[l-(l+^) 2 /3 2 /4]' 



(7) 



For fi s 0.3, this expression can be approximated as A w 
AqB^ within 1% accuracy in the range -0.1 < < 0.7. Here, 
A) = 2/x(l - Us), B = 7(1 + 0), and ip = 2.7, 2.9, and 3.1 
for fi s = 0.3, 0.4, and 0.5, respectively. 

4. BULK VELOCITY IN A PAIR-DOMINATED FLARE 

We now estimate the expected bulk velocity for a pair- 
dominated flare using a simple toy model. The power con- 
sumed by the e blob is the sum of two parts: the primary 
injected power due to dissipation of magnetic energy, Ldiss ~ 
L - L s , and the reflected radiation intercepted by the blob, 
f=s L s . Ldiss is probably injected in the form of relativistic par- 
ticles which are supposed to share their energy and momentum 
with the thermal e plasma, most likely due to collective effects 
or due to synchrotron self-absorption (Ghisellini, Guilbert, & 
Svensson 1988). Let $ d iss be the total vertical flux of the in- 
jected energy. The injected particles may accelerate or deceler- 
ate the thermal plasma depending on the angular distribution of 
injection. The reflected radiation always tends to accelerate the 
coronal plasma away from the disk. 

The energy streaming into the blob per unit time equals 
Ldiss + L s = L, and the corresponding net energy flux equals 
$diss + &s = The energy consumed during time dt in- 
creases the blob momentum by dp = (&/c)dt and the iner- 
tial mass by dm = (L/c 2 )dt. The net acceleration is d0 = 
d(p/mc) = (<j> — 0L)dt/mc 2 . The e ± plasma accelerates on 
a short time-scale (see §2), and the bulk velocity should relax 
to an equilibrium value for which the net acceleration vanishes. 
The equilibrium velocity is determined by the equation 



$-0L = O. 



(8) 



This equation also expresses the condition that the total energy 
flux vanishes in the comoving frame, $ c = 7 2 (<1 ) — 0L) = 0. 

We do not possess a detailed model of injection in the blob 
and do not know the feedback of the plasma bulk velocity on 
the angular distribution of the injection. Therefore, we consider 
two particular cases: 

(a) Ldiss is isotropic in the lab frame 



1 +Ms 



Ls 
L 



(9) 



On the other hand, L s /L is determined by equation (5). Equat- 
ing (5) and (9), we get an equation for the self-consistent bulk 
velocity and find ~ 0.1 for < ^ s ^ 0.5. 

(b) Ldiss is isotropic in the comoving frame 

In this case, the energy flux associated with Ldi SS vanishes in 
the comoving frame and the equilibrium velocity is determined 
by the soft radiation only, = $ S /L S . This yields = 1/2 for 

fi s = 0, and = 3/4 for fj, s = 1/2. 

A proton fraction exceeding m e /m p would increase the 
plasma inertia, and then the bulk velocity may fall below the 
equilibrium value. Pairs will tend to stream through the heavy 
proton component and a plasma instability may be initiated. In 
the stationary case, the stream velocity establishes itself just at 
the instability threshold which may be expected to be compara- 
ble with the thermal electron velocity, <~ c/2. 

5. CONCLUSIONS AND DISCUSSION 

The hard state of accreting black holes may be explained as 
a state in which a large fraction of the luminosity is released in 
a magnetic corona atop a cold accretion disk. The magnetic en- 
ergy is probably generated by the magneto-rotational instability 
in the disk and then transported to the corona due to buoyancy. 
The soft state may be observed when the coronal activity is sup- 
pressed and most of the energy is dissipated inside the disk. 

The energy release in the corona is likely to proceed in com- 
pact bright flares, where the local radiation flux strongly ex- 
ceeds the average surface flux from the disk. If the flaring 
plasma is comprised of pairs, it should be accelerated away 
from the disk by the pressure of the reflected radiation. The 
bulk velocity is then expected to be in the range 0.1 ^ ^ 0.7. 

The radiative acceleration of a pair plasma is one possible 
reason for bulk motion in the flares. Even if the flare is dom- 
inated by a normal proton plasma, the heating may be accom- 
panied by pumping a net momentum into the hot plasma at a 
rate <~ L /c. The transferred momentum per particle per light- 
crossing time, rb /c, is ~ lm e c. As the flare duration to 3> ry,/ c, 
the protons may acquire a momentum comparable to m p c. Note 
that magnetic flares where hot plasma is ejected towards the 
disk are possible. This case is formally described by equations 
(1-7) with < 0. 

The impact of the source velocity on the observed reflection, 
R (eq. [3]), and the Compton amplification factor, A (eq. [7]), 
is summarized in Figure 1 . Comptonization in the source pro- 
duces a power-law X-ray spectrum. The photon spectral index, 
T, is related to A by 



r w 2.33(>1 - 1)" 



(10) 



where 5 w 1/6 for GBHs and S w 1/10 for AGNs. This simple 
formula approximates within a few percent the results of calcu- 
lations we performed by using the code of Coppi (1992) (see 
Beloborodov 1999 for details). Combined with equation (7), 
this yields the dependence T(0). Assuming a typical [i s ~ 0.5, 
one can approximate T(0) as T w 1.9L? -0 5 for GBHs and 
T « 2L?-°- 3 for AGNs, where B = 7(1 + 0). 
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From the amount of reflection in Cyg X-l, R ~ 0.3, we 
infer a bulk velocity of the emitting plasma (3 w 0.3 for a 
system inclination <~ 30° (see the bottom panel in Figure 1). 
On the other hand, from the reported intrinsic slope r <~ 1.6 
and the corresponding A <~ 10, we again infer the same 
(3 w 0.3. Thus, (3 ~ 0.3 resolves the two problems of the 
disk-corona model mentioned in the Introduction. The beam- 
ing of the corona luminosity, Lx, also leads to the weakness of 
observed reprocessed blackbody component in the spectrum, 
Lbb ~ xRLx ~ Lx /4. Note that the inferred bulk veloc- 
ity is comparable with the thermal electron velocity, ~ c/2 for 
fcT-lOOkeV. 

The ejection model may apply to both GBHs and AGNs. The 
velocity of ejected hot plasma, (3, may vary (and in some cases 
it might be (3 < 0). An increase in (3 leads to decreasing R and 
T. A correlation between R and T is observed in the spectra 
of GBHs and AGNs (Zdziarski, Lubihski, & Smith 1999) and 
it is well reproduced by our model (Beloborodov 1999). Note 
also that fast outflows may manifest themselves in optical po- 
larimetric observations of AGNs (Beloborodov 1998b). 
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FIG. 1. — Top panel: The Compton amplification, A, of the reprocessed 
radiation as a function of the plasma bulk velocity, /3, (see eq. [7]). A repro- 
cessed fraction \ = 0.85 is assumed. The parameter p s describes the flare 
geometry (see eq. [4]). Bottom panel: The apparent amount of reflection, R, 
for the disk inclinations 9 = 0, 60°, and 90° (see eq. [3]). The observed R 
may be further reduced as a result of the scattering of the reflected radiation in 
the blob. The dotted line shows the expected A and R for (3 = 0.3. 
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